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ABSTRACT 


The results of three quasi-objective techniques for 
enhancing the accuracy of pen out cloud cover forecasts 
in sparse data areas,using satellite data and numerical- 
ly-produced analyses and prognoses, are presented. All 
three of the techniques involve the modification of the 
analyzed 500-mb relative vorticity, as based on the sat- 
ellite-observed cloud patterns and the associated 
changes to the prognostic fields of 500-mb relative vor- 
ticity. The final 24-hour cloud-cover forecast is based 
on the manually-modified relative vorticity pattern and 
three different advecting currents. These advecting 
currents are based on the prognostic 500-mb height field, 
modified and unmodified, and the prognostic 500-mb space 
mean field. The study was conducted over a large area 
of the North Pacific Ocean for the period 9 March to 


SP Apri? 1972. 
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2 ae INTRODUCTION 


The objective of this study was to investigate the 
possibility of increasing the accuracy of existing cloud 
cover forecasting techniques in sparse-data areas by us- 
ing satellite data in conjunction with numerical analyses 
and prognoses. In particular, an improvement was sought 
on Bittner's [1] quasi-objective method of using the ad- 
vection of relative vorticity on the 500-mb surface to 
determine the ee eee of cloud cover. 

Bittner's quasi-objective method uses current satel- 
lite-observed cloud patterns to modify the analyzed pat- 
terns of relative vorticity. The difference between the 
initial and modified relative vorticity patterns is then 
subjectively applied to the 24- or 36-hour prognostic 
patterns of relative vorticity. The primary areas of 
cloud cover are then related to the areas of positive 
relative vorticity advection (PVA) as determined by the 
unmodified 24- or 36-hour prosnostic 500-mb height 
fields (Method 1). 

Theoretically, an improvement in the cloud cover fore- 
casts could be made if the prognostic 500-mb height 
fields were modified by the amounts associated with the 
changes in the prognostic relative vorticity fields prior 
to determining the areas of positive relative vorticity 


advection (Method 2). 








milso, the theoretical development of Fjortoft [2] in 
manually producing a barotropic 24-hour forecast of the 
500-mb height values indicates that using prognostic 
space-mean fields to determine the areas of positive 
relative vorticity advection might show an improvement 
in the verification of 24-hour cloud cover forecasts 
(Method 3). 

Therefore, this study was specifically undertaken to 
determine which of the above three methods are the most 
accurate for forecasting 24-hour synoptic-scaled ee ua 
patterns. All three techniques were used simultaneously 
to produce cloud-cover forecasts over the North Pacific 
Ocean from 20°N to 50°N and from 170°E to 125°W. The 
forecast period extended Peom Fovaren to Lor April 1972. 

Modifying the 500-mb height prognoses in addition to 
the relative vorticity prognoses took slightly longer 
than the other two methods, but the resultant improvement 
in the forecasts seemed to indicate that the small amount 


of additional work was well worth the effort involved. 





If. BACKGROUND 


The oceanic areas of the world are relatively void 
of conventional meteorological observations, especially 
above the earth's surface. Over the middle-latitude 
areas the synoptic-scaled cloud patterns attain their 
best development. Fortunately, these cloud patterns, 
which are so revealing of synoptic systems to trained 
meteorologists, are easily observed with satellites. 
Consequently, for the past several years, meteorologists 
have striven to incorporate the vast quantity of synoptic- 
scaled information obtained from satellites into useful 
models for improved forecasts. Further, there has been 
considerable effort to incorporate this satellite-derived 
data into numerical analysis schemes, as Nagle and Hayden[7], 
Nagle and Clark [6!, Mantei and Workman [5], Zeigler (8] 
and Glaes [3]. Recently Bittner [1] has illustrated a 
method of using these satellite-observed synoptic-scaled 
cloud patterns to manually improve the numerical analysis 
of the 500-mb relative vorticity field. He then carried 
this improvement forward 24 to 36 hours by manually adjust- 
ing the numerically-produced prognostic 500-mb relative 
vorticity field. Finally, he converted this modified prog- 
nostic fieid into cloud and precipitation forecasts. 

In an equivalent barotropic numerical model of the 


atmosphere it may be easily shown that the sign of the 
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advection of relative vorticity at the level of nondiver- 
gence (LND) is the same as the sign of the vertical motion 
through this level. In fact, in this simple model of the 
RengeeNere the sign of the vertical motion is the same 
throughout a given column of air, if terrain effects are 
neglected. Therefore, positive advection of relative vor- 
ticity at the LND produces positive (upward) vertical 
motion throughout the column of air. In the real atmosphere 
the sign of the advection of relative vorticity, or geo- 
strophic relative vorticity, on the 500-mb pressure surface 
is highly correlated with the sign of the synoptic-scaled 
vertical motion through this surface. That is, if the 
advection of geostrophic relative vorticity on the 500-mb 
surface is negative then the synoptic-scaled vertical motion 
in the vicinity of the 500-mb level is also negative (i.e. 
subsidence) almost all of the time. 

If hygroscopic nuclei and moisture are present within 
the atmosphere then the formation and dissipation of 
clouds are largely determined by the sign of the vertical 
motion. That is, with the above conditions, positive 
vertical motion on a synoptic scale will produce synoptic- 
scaled cloud patterns when the time scale is sufficient 
to cover the phase lag between initial lifting and subse- 
quent condensation. With a given moisture distribution, 
the stronger the advection of vorticity the shorter the 
time lag between the vertical ascent and its tell-tale 


signature, the synoptic-scaled cloud patterns. 
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For a given analysis or prognosis of the height field 
of a constant pressure surface the relative vorticity may 
be obtained diagnostically from an equation which relates 
the mass field to tne wind field. The exact form of this 
diagnostic equation depends upon the assumptions made 
concerning the relationship between the wind field and 
the mass field. Fjortoft [2] demonstrated that the height 
field of a constant pressure surface may be smoothed by 
using a space differential operator called the Laplacian 
or del-squared operator. This smoothed height field is 
more commonly known in meteorology as the space mean. The 
Space mean represents the initial height field after the 
shorter wave lengths have been selectively removed. The 
exact short wave length which is completely removed from 
this field depends upon the grid length selected for the 
finite difference equation used in the approximation of 
the analytical equation. Fjortoft proved that, except for 
a latitude-dependent weighting function, the geostrophic 
relative vorticity field is obtained when the initial 
one field is subtracted from this space-mean height 
field. 

The height field of a constant pressure surface which 
is obtained through numerical analysis or DEE BRO OTE may 
be subjected to the scale separation technique of Holl [4]. 
In effect the initial height field is repeatedly smoothed 


with a Laplacian operator until a residual field called 
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the SR field remains. This SR field may be roughly con- 
strued to be equivalent to the space-mean field described 
above. When this SR field is subtracted from the initial 
height only a disturbance field, called an SD field, re- 
mains. The SD field roughly represents the combined effect 
of all the numerically representable short waves with wave 
lengths shorter than wave number six. Likewise the SR pat- 
tern approximately depicts wave numbers zero (polar vortex) 
through six. Note that the field which approximates the 
geostrophic relative per et et ty is the space-mean minus the 
initial height field while the SD field is minus the SR 
(space mean) plus the initial height field. That is, when 
the pattern of the SD analysis is used to approximate the 
pattern of the geostrophic relative vorticity, the sign of 
the SD centers must be reversed. See Appendix A. 

Consequently, the above concepts led to the development 
of quasi-objective rules for relating synoptic~scale clouds 
to the 500-mb SD field. Bittner [1] outlined these and for 
convenience they are listed in Appendix B. Bittner's 
method consists of using the most recent satellite-observed 
cloud patterns to modify the current 500-mb SD analysis. 
The resulting difference field between the initial and modi- 
fied SD patterns is applied to the prognostic SD field to 
produce modified forecast SD patterns. 

The particular step of modifying the SD patterns of the 


prognostic fieid based on the modifications applied to the 
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Current SD patterns appears to be the most subjective of 
all the procedural steps of Bittner's method. If there 

is no energy transformation occurring between the short 
waves represented by the SD field and the long waves por- 
trayed by the SR field during the period of the prognosis, 
then the modifications of the prognostic SD patterns are 
rather straight forward. However, for a given cyclonic 
Semareicyclonic system within an initial SD field, the 
greater the forecast change of development of a short-wave 
feature the less reliable the modification which is applied 
to the SD system of the prognostic field. That is, the 
more an SD center changes during a prognostic period, the 
more difficult it becomes to propagate an initial modifi- 
cation of the SD center to its prognostic counterpart. 

The modified SD prognosis is then used in conjunction 
with the unchanged 500-mb height prognosis to derive a 
subjective cloud and precipitation forecast valid at the 
time of the numerical height forecast. Bittner makes the 
following Pinpiioce relative to his synoptic-scaled 
cloud pattern forecasts: 

ih) Petes. cate vertical motion fields are highly cor- 
meiamedewith the characteristic distribution of extratropi- 
cal, synoptic-scale cloud patterns as seen by satellites. 

2) In turn, the vertical motion patterns are closely 
related to the advection of relative vorticity at 500 mb. 

3) An SD field is equivalent to a relative vorticity 


field. See Appendix A. 


14 






4) Any errors in the SD field based on the initial 
500-mb analysis are propagated during the numerical prog- 
nosis to produce similar errors in the SD field based on 
the prognostic 500-mb heights. Hence the error field is 


assumed to be conserved during the prognosis. 
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III. DATA SOURCES AND SELECTION OF STUDY AREA AND PERIOD 


Momarea Of the North Pacific Ocean extending from 20°N 
to 50°N and from 170°E to 125°W was selected for investiga- 
tion for several reasons. 

This area was selected primarily because it is both 
analyzed numerically by FNWC and has almost continuous day- 
light picture coverage by the ATS-1l1 satellite. ATS-1l is 
the designator for the geostationary Applications Technology 
Satellite positioned near the equator of the Pacific Ocean 
since November 1967. Its spin-scan cloud camera takes time- 
lapse photographs approximately every 25 minutes. The re- 
sulting pictures cover a region from about 50°S to 50°N over 
a span of 100 degrees of longitude centered at 150°W, with 
a resolution of two miles at the center. 

The subject ATS-l photographs were generally available 
locally from FNWC on a real-time basis Monday through Friday. 
Since the Saturday photograph was needed for verification 
of meey's forecast, the study period was usually limited 
to four days per week, namely Monday through Thursday. Some 
of these photographs and others that were missing on occa- 
sion were supplied by the Naval Air Systems Command's Pro- 
ject FAMos.? The FNWC computer products were also available 


locally on a real-time basis. 





1 ; ee 
Presently part of the Environmental Prediction Research 
Poicterty, Menterey, California. 
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The period selected (March and April 1972) not only af- 
forded a real-time basis for investigation but it is also 
the time of year when large, well-defined systems charac- 
terize the North Pacific Ocean. These systems readily lend 
themselves to the modification techniques evaluated during 


this study. 
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IV. PROCEDURES 


The research objectives of the procedures described 
memliow were to extend Bittner'’s work [1] and to outline an 
approach which would be practical for Fleet use. At the 
time when the FNWC 0000 GMT numerical analyses are received 
by the Weather Centrals and Facilities, the prognostic 
fields are not yet available. Therefore, to expedite 24- 
hour cloud forecasts from 0000 GMT it would be best to use 
the 36-hour prognostic height fields, as based on the pre- 
vious 12-hour analysis, and the current 0000 GMT analysis. 
Bittner indicated his approach is good for 24- and 36-hour 
forecasts. Since he effectively assumes that there is no 
energy transformation occurring between the SD and SR 
fields for 24 to 36 hours, then a change in an SD pattern 
may be propagated forward from any time step within this 
forecast period. That is, a change in the SD pattern at 
plus 12 hours from the time of the forecast may be carried 
forward to plus 36 hours. Therefore, in these procedures, 
changes noted in the 0000 GMT analysis of the SD patterns 
were incorporated in the prognostic SD patterns which 
verified 24 hours later, but which were made from the 
previous 1200 GMT analysis. 

Another point which may be mentioned in this regard is 
that within a sparse upper-air data area such as the North 


Pacific Ocean, the difference between the 12-hour prognosis 
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of the 500-mb level and the verifying 500-mb analysis may 
be quite small. 

The procedures in this study deal almost exclusively 
with FNWC's scale separation products reproduced on a 
1:60,000,000 polar stereographic projection. This required 
one scale change--that of manually transferring the cloud 
features of the ATS-1 photograph to this scale. FNWC is 
currently testing a program that will convert their pro- 
ducts to the scale and projection of the ATS-1l1 photograph. 
When this is perfected, all future studies can be carried 
out On the ATS-1l scale to eliminate any errors in the 
Manual transfer process. 

In general, this study consisted of determining a modi- 
fied SD field and projecting these modifications on the 
prognostic SD field. Then by using three advecting fields 
(FNWC's Z and SR PROGS, and a hand-modified FNWC Z PROG, 
called Z'PROG), three cloud cover forecasts were determined 
for verification. The detailed procedures are as follows: 

1) Each day the ATS-1 photograph closest in time to 
0000 CMT WVIsmsewecleaq = Tor Wse. “Mie™=major cloud Treatunes 
were traced onto a matching grid (same scale and subpoint 
as the photograph), carefully delineating as covered those 
areas with greater than 50% cloud cover and as open those 
areas with less than 50% coverage. 

2) The cloud patterns were then transferred by hand to 


@ polar-stereographic projection of 1:60,000,000 scale. 
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mare was a difficult and tedious procedure at first but 
after some practice, the time was decreased to less than 
15 minutes per map. 

3) FNWC's SD and SR analyses were put in register 
with the transferred cloud patterns. In areas of major 
cloud features such as fronts, cyclones, etc., the SD 
patterns and isolines were modified in accordance with 
the rules specified by Bittner [1]. These rules are listed 
in Appendix B for convenience to the reader. Further, areas 
of greater than 50% cloud cover were assumed to be associa- 
ted with positive vertical motion as indicated by positive 
vorticity advection using analyzed SR winds as the advecting 
field. The SD lines were redrawn to indicate the required 
advection in order to support the observed cloud patterns. 
Since the zero lines are not produced on the SD charts, 
these lines were drawn by hand. The result of the above wis 
a modified SD chart, hereafter referred to as the SD' 
analysis. The SR fields were not modified. 

4) Subjectively, the same adjustments that were made 
to ire SD analysis were then applied to the 36-hour SD 
prognosis, hereafter designated as the SD' prognosis. Since 
the location of zero-valued SD isolines are somewhat arbi- 
trary, care must be used to surround the same centers with 
the same zero line on both the analysis and the prognosis. 

9) The SD’ prognosis was then graphically added to 
the SR prognosis to produce a modified 500-mb height prog- 


nosis, hereafter indicated as a Z' prognosis. 
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6) The Z' prognosis was put in register with the SD' 
prognosis and used as the advecting field. Clouds were 
forecast in all areas of negative SD advection (positive 
Pebativesevorticity advection). Clear skies (less than 507% 
covered) were forecast for areas of positive SD advection 
(negative relative vorticity advection). Note that clouds 
(greater than 50% covered) were assumed to be located only 
in areas of positive relative vorticity advection in this 
Memdy. This is a slight deviation from Bittner's forecast- 
ing procedures for clouds as noted below. 

7) To produce the second cloud cover forecast, steps 
1) through 4) are the same but step 5)is omitted since the 
SR prognosis was used as the advection field for the pur- 
pose of step 6). 

8) To construct the third cloud cover forecast, steps 
l1)through 4) are the same but step 5) is omitted since the 
unmodified Z prognosis was used as the advecting field 
(Bittner'’s approach [1]) for the purpose of step 6). 

9) Steps 1) and 2) above were then applied to the 
next oe ATS~1 photograph (current plus 24 hours) to 
produce a verifying nephanalysis. 

10) The verifying nephanalysis was separately put in 
register with the three forecasts. Each of the 351 grid 
points for each of the three forecasts was checked and 
compared with the corresponding grid point of the verify- 
ing nephanalysis. Grid points which were observed clear 


but forecast cloudy were denoted by a "2". Those which 
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were both forecast and observed cloudy were indicated with 
a "O". Next, the points which were observed cloudy but 
forecast clear were denoted by a "1". Finally the grid 
points which were not covered by clouds on the verifying 
nephanalysis nor forecast to be cloudy were represented by 
a "3", At this point the technique which gave the best 
results could be readiiy determined by totaling the number 
Me correct forecasts ("0" and "3") for each forecast day. 
Tables I, II and III give the raw scores for the advecting 
fields of the Z' PROG, SR PROG and Z PROG, respectively. 
The sixth column represents the number of correct forecasts 
while the seventh column denotes the number of incorrect 
forecasts. 

11) Finally, the National Meteorological Center's 
threat score was computed for each of the forecasts. a 
explanation of the threat score is presented in Section VI. 

As indicated earlier, the manner in which the clouds 
were determined by this author is not strictly the same as 
that used by Bittner [1]. In addition, he applied the rules 
of Appendix B in reverse to indicate some geometric detail 
in the cloud patterns, as spiral clouds, vortices, etc., 
whenever they were deemed necessary as determined by the 
SD" prognosis. He also used various eae ines eee rules 
as a guide for forecast cloud tops, areas and amounts of 
precipitation, etc. Since this study was conducted to try 
and improve Bittner's basic premise, these refinements were 


not included. 


22 








9€T Gaz, adelaAy 








902 G7T L9 VET SL ai 7L/ Ye 
Cor 96T 7 1 COT Ca aS Ct] @ 
6ST Cin LZ €6 COT 99 CH 
99T VST €% TOT Tear C9 TT/% 
E7T 87Z OS 99 SLT LS Onley 37 
OOT TSie mA T9 Car 6€ 6/% 
€OT 8772 T4 cy EO? T9 8/% 
6L rd! Ce T? OSZ Be L/% 
TT OTZ Ga 86 88T Cv 9/47 
se 1 Caz TZ VOT Cec 77 c/% 
ZUT 602 T? TOT 89T TY 7/7 
TOT OS % T8 79 69T 6€ Lay 
C€T STZ L9 C8 TST OS C7/E€ 
Lom vTZ 6, ae Z9OT cS 7Z/€ 
UW 76, JE (sec YT T 79 CaET 09 Ciay © 
Cie 6LT C6 9TT 78 9¢ CORE 
947T Gc L6 76 80T ZS LAG 
Gor 972 96 pL O€T TS 9T/E 
Car 6LT 8S re 1 Hear Ov CT/E€ 
FE, 1 777 60T 8% CET 6S OT/E 
sles ees, LWD 00 
a ="'ly, ue + On wf. SEO LSOd LON uOu SO S40 LON 7 GT. 
= SLSO4 = SLSO4 aQNV LSOd LAG Sao QNV LSsod Lid Loe our L/eIed 
VIVIAOOUL VDIAAOY :XLBIATI >SpnoTg >SpnoTg :SpnoTy) UOFABIYVTI FAIA 





(‘sqjufod pyt23 yo xaqunu aie sontTea iazeptnaey) “PTOTFJ Sup Aooape 
9u2 SB 9NOUd,Z Fy? Butsn sAsesaxroF pnotTos jZufod-ptias fo uofAeoFFFI9A “Il VIGVLI 


23 




















c71 SG aSZeIsAYy 
Outed E77 89 VET EL 92 te 
"COT L6T 61 eae SIT as ET/% 
Gar 002 8Z 88 Ger €9 Cle 

TST OLE cE COT CET 91 ioe 
OG EG 0% 61 T9T ie OT/% 
TOT O52 Ly 09 CAT. TY? 6/4 
SIT E€Z cE LY 86T TZ 8/4 
76 L Sie OZ OS LEZ A] L/% 
99OT L8T Lee Zou 09T ZY 9/4 
9ST C6T Tec Gor oT 7Z G/¥% 
oat 002 €% VTT EGA Le 9/4 
8TT EEZ 88 Ez C7T mh, T/4 
tet On OZ 90T Oe T CE GZ/E 
6ST Z6T ws 90T O7T GG VZ/E 
oon 672 61T €9 Ore 6S EZT/E 
69T 78ST 96 Cit 98 AS CONE 
Eve OTZ 96 C6 ae 9% LT/E€ 
cae 972 88 89 SET LS OWE fae 
98T COT Gc O71 OTT va) CT/E 
Caer 6E72 eat as VET 09 OT/E 
at al pp LAD OO 
| ae ae ee oe wlu SO LSOI LON 10; SSO Sa0 LON 7L6T 

=~ SLSO4 =e SLSOd aGNVY LSOd Lng Sao aNV LSoOd LN@ LSod aut L/37ed 
FVIAAOOUT 4990LL09 >LCITO >SpnoTo >SPpnoTtgo >SPNOTO UOTABPOTFTFAVA 


2u2 SB 9OUd US FuA Butsn syAseosaxoF PpnotTo yupfod-pfas fo uvoFARBdI;FTAIA 


(*squtrod p}yx3 fo xtoqunu o1e SanTeaA zeTNGe_L) 





‘“PLTeTFJ Sut Ao.oaape 
“Di gagva. 


24 








Bani 702 . aZeIBAY 











C07 87T OZ 97T Sl al yT/% 
89T CST SZ Git cOT ES ET/4 
OST LOZ CE 76 99T 8S 7T/ 4 
7ST LOT O04 OLE Lot va hr / 2 
99T C8T LE 06 84 T 92 OT/% 
Sie 9€72 GL Tez. 79OT vay] 6/4 
en T 077 LS 99 EST C4 8/4 
66 ~ Siz ZZ cS 0€Z yy L£L/% 
94T COZ 6Z ee 9ZLT Ge 9/4 
Le TE Ode EZ 90T L6T CZ G/% 
oS. 66T ve OTT SOT 74 b/ 
eo | OZaZ EL 78 ay T 64 T/% 
79T 1 Or 89 66 6€T Cy G7/E 
mS 1 002 GS 66 c47T ZS wcole 
87T €7Z 9TT 89 LOT 09 SAAS 
Gt 9LT GOT 02 t aA cS GCs 
LET "TZ Gol L6 at T 04 Bl/¢ 
671 Zoe 76 08 O€T 64 9T/E 
O6T TOT ES Cv7T 80T Cv CT/E 
Tea 7 oe E OT C9 wer 79 Ot /e 
ai Ps ; ii (ar LWO 00 
nC ie, a ee ncnesdO. soa LON Oi Sao SaO LON Coot 
= SLSO4 =) S G04 GNV LSOa LN@ Sao aQNV LSOd LA@ LSoOd sUuTy/a ied 
RVIAANOOUT QVVBILTAAIND >A ROTO >Spnotg) > Spnoto >SPpnoto UOFAPIOTFIFAVA 
(*squyzod py23 Fo Aaqunu aze sanTeaA AzeTNge_L) "pTeTy sutjyooape 


29u2 S@® 9NOUd Z FU Sufsn syAseosxAoJ pnoTO Auyod-pyzaz fo woxFIReocFsFtAoA “III AMAVI 


2D 





Vqe soho ruD. OF MODIFICATIONS TO ANALYSES 


Piemsynoptuestime of /0000 GMT 10 April 1972 amd its 
associated analyses and prognoses were randomly selected 
as an example of the procedures presented in the preced- 
ing section. Figures 1] through 18 graphically illustrate 
the following discussion. 

Figure l is the ATS-l photograph of 2343 GMT 9 April 
1972 which was selected to be representative of the cloud 
cover for 0000 GMT 10 April 1972. The major cloud fea~ 
mumes (areas of > 502 cloud cover) are outlined for 
Clarity. 

In general, SD patterns moving into the western sector 
of the verification grid were modified very little because 
the cloud patterns are not well-defined near 170°E in the 
ATS-1l pictures. This sometimes led to a large error in 
the verification of cloud forecasts west of 180°. There 
was usually good definition of the cloud features and, 
therefore, pe cuee modifications of the SD patterns in 
the extreme eastern section of the grid. These eastern 
SD features and associated cloud features sometimes moved 
out of the area of the verification grid during the time- 
span of the prognosis. Some of the centers of the SD 
analyses in this near North American coastal area were 
extensively modified. For those systems which did not 


appear within the verification area at the verifying time 
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of the prognosis, the overall goodness of the associated 
modifications and the resultant cloud forecasts are un- 
known. 

Figure 2 shows a nephanalysis (areas of > 50% cloud 
cover as derived from Figure 1, on a grid of the same 
scale. 

Figure 3 is the nephanalysis after it has been trans- 
ferred to a polar stereographic projection with the SD 
and SR fields superimposed. The SD pattern is illustrated 
with dashed lines and the SR field with solid lines. Note 
the strong SD low center at 43°N and 135°W. Such a low 
value for the SD center suggests a cyclone at 500 mb. 
Since this center is embedded in an SR trough, as indicated 
by the solid lines, there is high probability of a closed 
center at 500 mb. However, the cloud patterns of Figure l 
indicate only slight cyclonic circulation with convective 
clouds in the area. This suggests decreasing the positive 
vorticity advection (negative SD advection) considerably. 

Figure 4 amet the new or modified SD field. Note 
the drastic changes made to the SD zero line. In the lower 
left portion of the figure, the zero line was relocated to 
the poleward edge of the frontal cloud band. (Rule l.a 
Appendix B). In general, SD troughs and ridges were re- 
drawn to coincide with major cloud boundaries. Note the 
considerable filling of the SD low center at 43°N, 135°W 
to decrease the negative SD advection. The center at 33°N 


180° was moved slightly northward to maintain the correct 


Peas 





perspective with the circulation center on the photograph. 
Meee tule 2.e Appendix B.) 

Figure 5 is the modified SD analysis of Figure 4 
(dashed lines) superimposed on the original SD field of 
Figure 3 (solid lines). The amount and direction of move- 
ment of all centers and isolines were noted and used in 
the modification of the SD prognosis. 

Figure 6 represents the 36-hour SD prognosis and its 
Mediiication (SD PROG). As such, it represents the 
changes made to the SD analysis 24 hours earlier, as shown 
in Figure 5. As can be seen from the figure, this modifi- 
cation is somewhat subjective. 

Figure 7 shows the graphical addition of the modified 
BD prognosis (SD" PROG) to the SR prognosis in order to 
produce a new height prognosis (Z’' PROG). This is the 
only additional step required over Bittner's method. 

Figure 8 is the Z' boentes te of Figure 7 superimposed 
on mire SD’ prognosis of Figure 6. Since the negative SD 
values Potrespoua te positive values of vorticity, as 
shown in Appendix A, the areas of positive vorticity advec- 
tion (PVA) : readily apparent. These areas are carefully 
delineated and are represented in the figure as clouds. 
Note that major cloud boundaries are determined by the SD 
troughs and ridges as defined earlier relative to vorticity 
advection and as stressed by Bittner. As mentioned in the 
previous section, Bittner inserted cloud-pattern detail 


(i.e. cloud vortices, developing cyclones, etc.) at this 
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Point. these modifications could be done by a knowledge- 
able synoptic meteorologist resulting in better GOmeianii ty 
of cloud patterns and increased accuracy of cloud fore- 
casts. 

Figures 9 and 10 are the same as Figure 8 except that 
the vorticity-advection areas are determined by using the 
SR unmodified Z prognoses as the advection fields, respec- 
tively. See procedures 7 and 8 of the previous section. 

Figure ll is the verifying ATS-1l photograph with areas 
fe? 202 cloud cover outlined. This picture also shows 
very little cyclonic motion in the northeast corner, justi- 
fying the decreased PVA of Figure 4. 

Figure 12 shows the boundaries of the > 504 cloud cover 
from the photograph of Figure ll. 

Figure 13 illustrates the verifying clouds as trans- 
ferred from Figure 12 and the PVA areas derived in Figure 
8 from the Z' prognosis. As mentioned earlier, the date 
for this case study was chosen at random; it just happens 
to be a below-average date. See Tables I, II and III for 
comparisons of raw values for 0000 GMT 11 April 1972. 

Figure 14 represents she Awetsags of the forecast, 

The 351 grid points of Figure 13 (every 2-1/2° latitude- 
longitude spacing) were examined and assigned to one of 
four categories. The grid points that were forecast to 
pe clotidyebut were observed to be clear ("2") totaled 65. 
Those that were both forecast and observed cloudy ("0") 


amounted to 141. Observed cloudy but forecast clear ("1") 
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were 101 points. And, finally, those grid points which 
Mere both forecast and observed to be clear ("3") totaled 
3. This amounted to 184 correct ("0 + 3") and 166 ("1 + 
2") incorrect forecasts, which yields a combined threat 
Bootes” of .67. See discussion in Section VI and scores 
in Table IV. 

Figures 15 and 16 show the verifying clouds superim- 
posed on the PVA area derived from the SR prognosis and 
the scoring grid, respectively. The number of grid points 
with correct forecasts ("0 + 3") amounted to 170 while the 
jm@errect total (‘1 + 2") was 18). The threat seore for 
this method was .59. 

In general, the scores for all three methods in this 
case study were below their respective averages. See 
Table IV for the averages. However, the relative differ- 
ences between the three methods are of the same order as 
the overall averages of each. That is the SR method of 
determining the PVA areas yielded approximately the same 
score as Been cee a iast 4 we (unmodified Z field) for 
determining PVA areas while the Z' (modified Z field) 
me«éhod gave Merce: than a 104 increase in the value of the 
combined threat score which was used to evaluate the cloud 


cover forecasts. 


mainte? score will be explained in detail in the following 
section. 
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VI. VERIFICATION OF RESULTS 


The National Meteorological Center's threat score sys- 
tem was used to measure the goodness of the forecasts. 
The score depends on the correct forecasting of cloudy 


(Ty )} and clear areas (Ty). They are calculated as follows: 


0 2 3 
‘oe aoe eae 
where: 

O = Those points both forecast and observed as 
eiowd.. 

1 = Those points forecast clear but observed as 
cloudy. 

2 = Those points forecast cloudy but observed as 
clear. 

3 = Those points both forecast and observed as 
clear. 


Using this approach, a perfect forecast of clouds re- 
Sults ina score of T 6 = 1.0 and a perfect forecast of 


clear skies results in a score of To = le OG a  COmlis a mid 1o 
these two perfect threat scores gives a total of 2.0 for 


an absolutely correct forecast over an area as follows: 


+ = ], 0 = 2, 
Le i iste) 2.0 


The results of the three forecasts prepared for each 
day are listed in Table IV. These scores were obtained 
using the above technique. A cloudy forecast score and 


a clear forecast score followed by a combined score are 
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peesented for each day of the study period. The overall 
averages of each method are listed at the bottom of the 
table. 

The results of Bittner's study of approximately the 
Same area of the North Pacific Ocean, but for different 
dates, are presented in Table V for comparison. The un- 
modified column represents threat scores obtained without 
the use of satellite pictures to modify either the initial 
or prognostic SD fields. The modified column represents 
the forecasts obtained after the initial SD fields have 
been modified based on the satellite observed cloud pat- 
terns and the prognostic SD fields have been modified 
based on the modifications applied to the initial SD 
patterns. This latter column represents a 15% increase 
in cloud-cover forecasting accuracy over the unmodified. 
column. Except for the additional refinement steps men- 
tioned earlier, this modified column is equivalent to the 
Z PROG column of Table IV. The additional refinements in 
forecasting cloud Sorex used by Bittner may account in 
large part for the relative accuracy differences between 
the two corresponding columns of Tables IV and V. How- 
ever, as noted earlier, these refinements may be easily 
added to the basic approach in determining the PVA-associa-~ 
ted cloud cover. 

The basic determination of the PVA areas by the three 


different methods used in this study are clearly comparable. 
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mreeresults in Table TY amdicate that the Z PROG and the 

SR PROG advection fields are about equal in their results. 
THe zZ" PROG advection fields indicate about a 10% improve- 
ment over the approach suggested by Bittner [1]. Thus it 
may be argued that if Bittner's approach is modified to 

mee sextent gof including the Z’ (viele the Z) pro@tfosis for 
the advecting field, before adding the required and obvious 
continuity of large-scale cloud patterns, the best 24-hour 


cloud-cover fields may be produced. 
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TABLE IV. 
PROG and Z PROG cloud-cover forecasting techniques. 





 compeardsoms of thaeatescemes of Z'PROG, 


Verification 

DATE /Time Z"PROG SR PROG Z PROG 

13) 7d rh i eo Ab + T, =T T + = T 
00 GMT © e 

3/10 295. 5 1 ae PG «48a 1.02 49 .46 .95 
39/15 neal. 2.5 . 66 A374 23 60 soto 622 Res 8 
3/16 <Srleee 43 94 wo. a 93 —0mr 442 92 
ay17 74 3oen 40 83 aa AO ~84 4e «=— 43 . 88 
31/22 Beis) Peis 69 3 256 .70 -30° «37 sO 
3/23 46 ~48 96 ey ~49 ~96 40 - 48 94 
3/24 , 5 Gm 82 eee 6 25 Hie. sSOm <27 iT 
3y7/25 mesnc)) fenesi 87 = 5 (Cress: aes oe. 32 A eal 
4/1 O72 245) wee £55, 43 98 <3) «36 ~ 89 
4/4 lee 63 ol 22 73 se .18 .70 
4/5 61 .14 rE 55) Seer 65 oor .15 ws 
4/6 5 Jee ATE oh e498 24 63 ome .16 7) ah 
4/7 «fous 22 98 a2. 8 90 709 .18 88 
4/8 ~67 .29 ~96 Ors 23 » 86 -62 ~ 34 ~96 
4/9 , 6 3 4 Cees. On 163. «43°. 8.06 -59 .38 97 
4/10 59 .29 88 275) ne cal 96 49—) 43 65 
4/11 046 21 .67 ensee 59 arlene. qales 59 
io -89 «150g 53 «6 69 ose 3 AY 
dyed HAE 1: 76 se ia Se sles <3iee «32 BHO 
4/14 Ay en? 5 26 Se Bylo 36" a2 59 
Average .56 31 87 pe ee 28) 80 -49 .29 IES: 
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Peer V. The results of Bittner’s study [1]. 

DATE UNMODIFIED MODITIED 

1969 Ty ti Ty = T T =F = T 
O07 3 woe, «46606 )| 1208 2 61 isis: 
10/9 <310 .63 093 oS -73.. Lo 
20/10 29 ~49 Aaviae: ot ~66 1.00 
10/13 24 ~58 Bree oe woe - 66 98 
10/14 n4 60 aad reas . 03 Asie) 
1037 1.5 m2 ~29 < oF soe soo. L..0:0 
10/16 2h saa oul mo, we 120s 
10/17 i220: <6eealesOS H26» 370 ~96 
10/18 Ze Oe wor ~so”6 Ch CU e O 
10/19 De naif See = ey oi a Oe 
10/20 ey) 7) Bs) ~ 89 eo 106 99 
deQy/ 2 1 ae wis oy, 5 68 +71 eS 
10/24 26 6.43 .69 [25 °6~=Cté«‘( CCS 793 
20/25 638. oe »99 57) (eo) mo 
10727 25 eo ~85 of 3 (se ere 2a 
10/30 26 ae wide: ~40 vos 1.03 
a7 4 7o2 oom ~ 89 Cee D / mo 
Average 29° fom 88 7o2e oo) 1 01 
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VII. CONCLUSIONS 


This study indicates there is merit in using the satel- 
lite-observed cloud patterns to modify the position and 
orientation of 500-mb relative vorticity patterns (the SD 
fields of FNWC) in order to improve 24-hour forecasts of 
cloud cover over oceanic areas, Betene: sig. in an earlier 
Study, came to the same conclusion. 

It appears that Bittner's basic approach can be improved 
by using both modified prognostic relative vorticity pat- 
terns and modified prognositic height fields at 500 mb. 
However, the refinements of cloud-pattern continuity de- 
scribed by Bittner should be included in the final 24-hour 


cloud pattern forecast for best results. 
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2343 GMT 09 April 1972 ATS-1 satellite photo- 
Area of > 50% cloud cover bounded 


Macure 1. 
gmaph (Subpoint: .2°S). 
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Nephanalysis of ATS-1l satellite photograph 


Bajoare 2), 
Area of > 502 


Gsirb petatee a)4 2343. GMT 09 April 1972. 
cloud cover bounded by @. 
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Figure 3. FNWC 500-mb SD (---) and SR (——) analyses, 0000 
GMT 10 April 1972, superimposed on the ATS-1 nephanalysis 
oe > 50% cloud cover (2343 GMT 09 April 1972). SD interval: 
30 m; zero SD line labeled. 





Figure 4. FNWC 500-mb SR (-—-——) and modified SD (---) 
analyses, 0000 GMT 10 April 1972, superimposed on the 
miSs-i Wephamalysis of > 5024 cloud cover (2343 GMI 09 

April 1972). SD interval: 30 m3; zero SD line labeled. 
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Piegure 5. 
analyses, 
SD line labeled. 


Figure Gi. 











500-mb FNWC SD (——) and modified FNWC SD (---) 
0000 GMT 10 April 1972: SD interval: 30 m; zero 
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36-hour 500-mb FNWC SD (-——) and modified FNWC 
SD (---) prognosis, verifying at 0000 GMT 11 April 1972: 
SD interval: 30 m3; zero SD line labeled. 
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Figure 7. Modified FNWC Z prognosis (-—-——) as derived from 
the FNWC SR (-*-*) and modified SD (---) prognosis, verify- 


Mme at OjM@O@P GMT 1) Apis loOWe: SD interval: 30 m; zero §D 
line labeled. 








Figure 8. Modified FNWC 500-mb Z ( 


noses with resulting areas of positive vorticity advection 


(€3 ), verifying at 0000 GMT 11 April 1972: SD interval: 
30 m; zero SD line labeled. 


) and SD (---) prog- 


40 











Figure 9. FNWC 500-mb SR (¢ ) and modified SD (---) prog- 
noses with resulting areas of positive vorticity advection 
(© ), verifying at 0000 GMT 11 April 1972: SD interval: 
30 m; zero SD line labeled. 





Figure 10. FNWC 500-mb Z (——) and modified SD (---) prog- 
noses with resulting areas of positive vorticity advection 
( C3 ) evry ime at 0000 GMT 21 April 1972: SD interval: 
30 m; zero SD line labeled. 
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Hoemre |i. 2s Gn 10 April POV 2esT>-l™earellite photo- 
graph (Subpoint: 2°S). Area of > 50% cover bounded by oP, 





Figure 12. Nephanalysis of ATS-1 satellite photograph 
(Gti podmeweemS), 2358 .GMT°LO “April .dO72y, Area of > 507 
cover bounded by @ . 
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Peeure 13. Nephana@lysis of > 30% clowd cover from ATS-1 
pPeotesraph (C2358°20r April 1972) ( GC and sreawof posi-- 


tive vorticity advection, (#55 ) derived from the modified 
36-hour FNWC 500-mb Z prognosis verifying at 0000 GMT 
feempril 1972. 
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Figure 14. Scoring grid for the cloud forecast based on the 
modified 36-hour FNWC 500-mb Z prognosis, verifying at 0000 
GMT ll April 1972. Legend: O-Clouds: forecast and observed- 
141 points; 1-Clouds: observed but not forecast-101 points; 


2~Clouds: forecast but not observed-66 points; 3-Clear: 


fore- 
cast and observed-43 points. 
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Figure 15. Nephanalysis of > 504 cloud cover from ATS-1 
photograph (2358 GMT 10 April 1972) (& ), and area of 
positive vorticity advection derived from 36-hour FNWC 
500-mb SR prognosis verifying at 0000 GMT 11 April 1972. 


Figure 16. Scoring grid for the cloud forecast based on the 
modified 36-hour FNWC 500-mb SR prognosis, verifying at 0000 
GMT 11 April 1972. Legend: O-Clouds: forecast and observed- 
135 points; 1l<~Clouds: observed but not forecast-105 points; « 
2-Clouds: forecast but not observed-76 points; 3~Clear: fore- 
cast and observed-35 points. 
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Figure 17. Nephanalysis of > 504 cloud cover from ATS-1 


photograph (2358 GMT 10 April 1972) (@ )}, and area of 
positive vorticity advection derived from the 36-hour 
FNWC 500-mb Z prognosis verifying at 0000 GMT 11 April 1972. 
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Figure 18. Scoring grid for the cloud forecast based on the 
modified 36-hour FNWC 500-mb Z prognosis, verifying at 0000 
GMT 11 April 1972. Legend: O-Clouds: forecast and observed- 
127 points; 1-Clouds: observed but not forecast-110 points; 
2-Clouds: forecast but not observed-74 points; 3~Clear: fore- 
cast and observed-40 points. 
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Figure 19. The relationship between cloud patterns and the 
500-mb SD field, as presented by Bittner [1] and illustrated 
by Zeigler [8]: A. Surface wave vs. SD field. B. Occluded 
system with 500-mb flow and SD field. C. Mature cyclone vs. 
oD fiem@deeD. Dissipating cyclome vs. SD field. E. -Polar jet 
axis and open-celled cumulus vs. zero SD contour. F. SD flow 
vs. cloud-associated PVA MAX. 
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APPENDIX A 
FNWC'S 509-MB SD FIELD AS AN APPROXIMATION TO 
THE RELATIVE VORTICITY FIELD 
The definitions and discussion below follow closely 
emeat in [1]. 


The 500-mb disturbance component is defined as 


. _ 2 2 
SDeog = 25997 S8sq9 = ~9S oa VO Z 6 9de (1) 
where: 
Ze00 = 500-mb height field 
SRe O09 = 500-mb residual (mostly wave numbers 1,2,3,4, 
5 and 6) 
SD. og = 500-mb shortwave field (mostly wave numbers 


7 and above) 


es 
I 


degree-of-smoothing parameter 
6S = numerical grid length 
The right hand side of this equation is a Laplacian of the 
Hemewt fletd over the Saouenunie interval a . 
Geostrophic relative vorticity is defined as: 
-~l 2 


C= gf Vi Zeng (2) 


where: g = acceleration due to gravity 


rh 
I 


Coriolis parameter 
By substituting Equation (2) into Equation (1), the follow- 


ing relationship is obtained: 


SD = -6S° fg f Td (3) 
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Equation (3) shows that the SD field is about the same as 
the geostrophic relative vorticity fieid (with a sign re- 
versal) over the smoothing interval used in the analysis. 
Thus, for all practical purposes, they can be considered 


analogous for both analyses and prognoses. 
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APPENDIX B 
SUMMARY OF RULES RELATING CLOUD 
PATTERNS TO 500-MB SD ANALYSES 
Following is a summary of rules relating cloud patterns 
to the 500-mb analyses [1]. 
1) The SD zero line should be drawn: (See Figure 19 for 

a graphical illustration of some of these rules.) 

a) Along the poleward edge of the cloud band associated 
with the polar front. 

b) Along or just equatorward of the cioud shadow mark- 

‘ing the polar jet axis; along or just inside the 
poleward limit of transverse lines marking the sub- 
tropical-jet axis. 

c) Along the line between two major vorticies which 
are connected by a continuous frontal band. 

2) An SD low center should be located: 

a) Five to seven degrees upstream from an open frontal 
wave, or at the rear edge of the head of a cloud- 
identified positive vorticity advection maximum 
(PVA MAX). 

b) Five degrees upstream from the center of a hook- 
shaped cloud pattern (rapidly developing cyclone). 

c) Within three degrees or less from the center of a 
complete spiral cloud pattern (mature cyclone). 

d) Directly over or slightly to the east of the center 


of an open vortex (dissipating cyclone). 
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a) 


>) 


6) 


e) 


An 


a) 


b) 


An 


a) 


b) 


c) 


The 


a) 


b) 


é.) 


SD 


a) 


b) 


Over the c@mter of a strong vorticity center in 


the cold air. 


SD trough should be drawn: 


Through the point where a frontal cloud band 
forks or fractures. 
Five to seven degrees upstream from a weak frontal 


wave (where a closed SD low is not justified). 


SD ridge should be drawn: 


In an occlusion, to coincide with the downstream 
major cloud boundary or the change from C (cumuli- 
form/ stratpirorm este MICO (stratitorm,cirriform). 
With anticyclonic curvature near the equatorward 
edge of a frontal cloud band. 

With @antacyclomic cunmvature over the peak of a 
frontal wave pattern. 

SD gradient should be maximized: 

Along the rear edge of a frontal cloud band. 

In the dry tongue of a hook-shaped vortex. 


At and poleward of the peak of a frontal wave. 


lines should be drawn: 


Parallel to lines of cumulus (CELLS OPEN) in the 
cold-air sector of a cyclonic disturbance. 
Parallel to cloud lines and bands in a dissipating 


vortex. 
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